On the reversal of quantum many body dynamics and its complexity 
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We demonstrate that via optimal control arbitrary time evolutions of many-body quantum systems 
can be reversed. The optimal reversed dynamics -contrary to standard time-reversal procedures- 
is extremely robust to external sources of noise. We suggest a relation between control complexity 
and the dimension of the manifold supporting the dynamics of a quantum system. We verify the 
validity of this conjecture in three different models elucidating the role played by integrability in 
this context. 

PACS numbers: 



In recent years, fast progress on the understanding of 
non-equilibrium dynamics of many-body quantum sys- 
tems has been spurred by unprecedented opportunities 
offered by cold atom quantum simulators as well as by 
the development of powerful numerical tools [l], [1| . In- 
deed, numerical methods have made it possible to in- 
vestigate the out-of-equilibrium dynamics of many-body 
quantum systems and to compare theoretical results with 
experimental data obtained in highly controlled and clean 
systems. Many interesting effects have been investigated 
so far Q including (just to give a few examples) quench 
dynamics , thermalization Q , quantum phase tran- 
sition dynamics B, effects of disorder ^ and of periodic 
perturbations [10, LUj both in fermionic and bosonic sys- 
tems [12]. These advances in manipulating and charac- 
terising non-equilibrium dynamics may result in the de- 
velopment of strategies to control the quantum evolution 
of many-body systems. This task, beyond the bounds of 
possibility until a few years ago, paves the way for the re- 
alisation of many-body state engineering where optimal 
control techniques [13] appear to be the ideal tool to use. 

Quantum optimal control has been applied only re- 
cently to quantum many-body systems, e.g. to adiabatic 
dynamics through a quantum critical point [14] or to the 
cooling of Luttinger liquids [l5|. An algorithm for opti- 
mal control particularly suited for many-body problems, 
the Chopped Random Basis (CRAB) optimization, has 
been recently developed and applied to state prepara- 
tion of strongly interacting cold atoms in optical lattices 
and spin systems [lf|. The theoretical study and ex- 
perimental implementation of optimal control strategies 
to quantum many-body systems poses in turn a number 
of important questions. While it has been shown how 
quantum optimal control can drive a system up to its 
quantum speed limit [I?], Eij •> ^ * s important to under- 
stand to which extent is it possible to control a quantum 
many-body system. Which resources are needed in terms 
of complexity, in particular in connection to the inte- 
gr able/chaotic nature of the system under investigation? 
And how efficient and robust will the resulting control 



Quenches 




FIG. 1: Dynamical scheme to optimally reverse the system 
dynamics: a system in the ground state is taken out of equilib- 
rium by multiple random quenches. Reversing the dynamics 
can be obtained in general via a time-inversion or by solving 
an optimal control problem. 



strategy be? 

In this Letter we address these questions by character- 
ising qualitatively and quantitatively the ability to reach, 
from a given input, a predetermined final state of a many- 
body system through the time dependence of some global 
couplings. We will in particular consider the problem of 
reversing complex quantum dynamics: a system is first 
taken from the ground state to a highly disordered state 
via random quenches (see Fig.^jK characterised by a high 
diagonal entropy Sd (see Ref. [19[), and then driven back 
to the initial state. We will first address the robust- 
ness problem, showing that while the obvious strategy 
of returning to the ground state by time-reversing the 
protocol is highly sensitive to external noise in the cou- 
plings [19| , an optimal protocol based on CRAB is in turn 
much more robust against errors, making it a good can- 
didate for experimental realizations. We then address 
the problem on a more general ground, elucidating for 
the first time the relationship between control complexity 
and dimension of the manifold in which the dynamics oc- 
curs. In particular, we demonstrate that the complexity 
is weakly dependent on the initial and final state prop- 
erties while it is strongly influenced by the integrability 
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FIG. 2: Diagonal entropy Sd for the Ising model (black line) 
and for the LMG model (red line) as a function of time during 
the disordering procedure (left side of the figure, logarithmic 
scale for the time) and during the optimization (right side of 
the figure, linear scale for the time) with N — 20, Tmax = 
100/ A, where A is the critical gap. Time is in units of J. 



of the system. 

These results bring together in a new perspective, op- 
timal control, complexity theory and thermodynamics, 
paving the way to further developments and investiga- 
tions. In particular, an interesting issue related to our 
study is the reversibility of closed many-body quantum 
systems dynamics, which might have intriguing conse- 
quences on a fundamental problem in physics, that is, 
the emergence of the arrow of time. Indeed, one can re- 
vert the dynamics by inverting the time propagator, as it 
is typically done in spin-echo experiments [2l| . This pro- 
cedure is however a highly non trivial task in a general 
many-body system and requires an enormous accuracy 
in the knowledge of the history of the dynamical process 
and control of the system: the minimum deviation from 
the exact path inversion has dramatic consequences [HI . 
In practice, very few systems are amenable to such opera- 
tions, since quantum systems are typically only partially 
tunable. The results presented here suggest that the in- 
trinsic complexity of inverting time evolution rests on the 
tensor-product structure of many-body quantum system 
Hamiltonians, that is, on the fact that in a non-integrable 
system the amount of information (or band- width) to be 
encoded in the optimal controlling field scales exponen- 
tially with the number of constituents. 

Optimal reversed dynamics- In the following we will 
focus on the following dynamical scheme (see Fig. [I]): i) 
we first initialize the system in an eigenstate (Sd = 0, sta- 
tionarity); ii) then apply a disordering quench process, 
generating high diagonal entropy Sd (out of equilibrium 
states); Hi) and finally steer the system back into the ini- 
tial state (Sd = 0) either using an optimized protocol or 
time-reversal. For our numerical computations we con- 
sider systems described by the following TV-spin quantum 



spin Hamiltonian: 

N N 

H = -J2 ~ r(t) J2 < ~ J * E w 

i,j i i 

For vanishing longitudinal field (J x = 0), this Hamilto- 
nian has two obvious integrable limits, the quantum Ising 
chain in transverse field, where = Jtf^+i, and the 
infinite-range quantum Ising model (or Lipkin-Meshkov- 
Glick (LMG) model [23]) when = J/N for i < j 
(Jij = otherwise). In the presence of a longitudinal 
field (J x 7^ 0), the quantum Ising chain loses its integra- 
bility 1241 . apart from at the critical point (in the scaling 
limit) [25j. From now on we set h = 1 and time is ex- 
pressed in units of J. 

Due to the various integrals of motion, the dynamics 
of the LMG model is effectively described in the sub- 
space generated by the Dicke states \S,S Z ), where S 
is the conserved total angular momentum and —S < 
S z < S are the allowed possible projections along the 
z-axis [26|. The ground state of the Hamiltonian be- 
longs to the subspace with S = N/2 and in the fol- 
lowing we are working within this subspace. We stress 
that in such subspace there are ~ N/2 Dicke states dy- 
namically accessible (corresponding to the allowed values 
of S z with the correct parity) and therefore the maxi- 
mal diagonal entropy scales logarithmically with the size, 
^ MG ~log(7V/2 + l). 

The Ising chain in transverse field for J x = can in turn 
be solved exactly through the Jordan- Wigner transfor- 
mation mapping the spins onto free fermions [27j . For the 
Ising chain the maximal diagonal entropy scales linearly 
with the size S T d ~ (N/2) log(2); indeed the dimension of 
the Hilbert space is 2 N I 2 due to parity conservation for 
the number of fermions obtained with the Jordan- Wigner 
transformation. 

For both models, we prepare the system in the ground 
state 1*0 (0)) = \GS) at large value of the driving field F 
- a fully polarized spin state along the positive z-axis. 
We then drive the system out of equilibrium performing 
a repeated quench between two values Ti and T2 of the 
control field T; each quench lasts a random waiting time 
Tmax ' r i-> where T max is the maximum allowed waiting 
time and Tj G [0, 1] is a uniformly distributed random 
number [19[. We verified that after a sufficiently large 
number of cycles the average Sd produced with this dis- 
ordering procedure is approximately independent of the 
amplitude |Ti — T2 1 and of the waiting time between two 
consecutive quenches. Finally we use optimal control to 
drive the system from the out of equilibrium state back 
to the initial state |^(0)), in a given time T to obtain 
the final state \ip(T)). The optimization is implemented 
through the CRAB technique [l6j: after making a guess 
for a possible return path To(t), we introduce a correc- 
tion of the form T(t) = To (£)/(£), where the function f(t) 
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FIG. 3: Infidelity I as a function of the intensity of the noise 
f in the Ising model for N = 20, with T(t) = T(l + f • r(t)), 
where F is the driving field in the absence of noise and r(t) is a 
random variable uniformly distributed in [—1, 1]. Black circles 
continuous (dashed) line: exact reverse dynamics T ~ 2 • 10 3 
(T ~ 90); red triangles continuous line: optimized reverse 
T-50. 



is expressed as a truncated Fourier series, i.e. 



m = i 



J2 k A k sin(v k t) + B k cos(v k t) 
A(t) 



(2) 



Here, fc = l,...,n/, ^ = 27rfc(l + r*fc)/T are random- 
ized Fourier harmonics, T is the total time evolution, 
rk G [0, 1] are random numbers with a fiat distribution, 
and X(t) is a normalization function to keep the initial 
and final control pulse values fixed. The optimization 
problem is then reformulated as the extremization of a 
multivariable function F(Ak, Bk,Vk), which can be nu- 
merically performed with a suitable method [l6| . 

In Fig. [2] typical results of this procedure are shown 
for the LMG model (red line) and the Ising model (black 
line): in the left half of the picture the disordering pro- 
cess is applied and the diagonal entropy increases reach- 
ing an average maximum value (S^ MG ~ log(7V/2) and 
S T d ~ 7V/21og(2)). In the right half of the picture the 
optimization phase is shown: even though the Hamilto- 
nian is only partially tunable (i.e. we do not allow for a 
sign reversal of all couplings) the control is able to steer 
the system towards the desired initial stationary state, 
reducing Sd to zero, without the need for information on 
the history of the disordering process. We verified that 
this holds for different system sizes and different total 
times T (data not shown) . 

As discussed in Ref. [19|, driving a system back to its 
initial state by a full time-reversal of a protocol is a pro- 
cedure extremely sensitive to very small noise perturba- 
tions. It is therefore natural to compare the effects of 
noise on a naively time-reversed protocol to those on an 
optimised CRAB protocol. In order to do so, we will con- 
sider a perturbed protocol T(t) = T(t)[l + r (£)£], where 
T(t) is the original protocol (including a sign change of 
the couplings for the time reversed one), r(t) is a random 
variable uniformly distributed in [—1,1], and £ is the in- 
tensity of the noise. In Fig. [3] we plot the final infidelity 




FIG. 4: Infidelity I as a function of the rescaled number of fre- 
quencies nf/N a in the Ising (upper panels) and LMG (lower 
panels) models, for the two transitions \C) — > \GS) (left pan- 
els) and \Ms) -> \GS) (right panels) with T = 50 (ising) and 
T = 100 (LMG) and N = 10, 20, 30, 40 (black circles, red 
squares, blue triangles, green triangles respectively). 



I as a function of £ for the exact reverse dynamics (black 
circles continuous line) and for the optimized reverse dy- 
namics (red triangle continuous line). The robustness of 
the optimized protocol emerges strikingly: a noise more 
than six orders of magnitude more intense needed to af- 
fect the optimal protocol to an extent comparable to the 
time-reversed one. In order to verify that the stabil- 
ity of the optimized protocol is not simply due to the 
reduced total evolution time, we repeated the analysis 
with a time-reversed protocol lasting a time comparable 
to that of the optimized process: the minimal time to 
obtain a maximal entropy state is around T ~ 90 (black 
circles dashed line), to be compared with the duration of 
the optimal process, T ~ 50. The results shown in Fig. [3] 
confirm that the optimal protocol is intrinsically more 
stable than the time-reversed one, making it an excellent 
candidate for experimental implementations. 

Control complexity- We demonstrated the feasibility of 
reversing a protocol connecting maximal entropy states 
with vanishing entropy states. We would like now to 
characterize the complexity of such a generic state-to- 
state conversion problem. Let us first give an opera- 
tional definition of complexity: within the CRAB algo- 
rithm we will measure the complexity in terms of the 
number rtf of Fourier components needed to solve the 
optimisation problem up to a certain target infidelity. 
Intuitively transformations from a state with maximal 
diagonal entropy (i.e. completely delocalized in phase 
space) to one with low diagonal entropy (i.e. well lo- 
calized) are expected to be more difficult than those be- 
tween localized states or between adiabatically connected 
states. It turns out however that the complexity of an op- 
timisation protocol depends only weakly on the choice of 
initial and final states. Let us illustrate this considering 



two different state-to-state transformations: from a max- 
imal entropy state to the ground state (\Ms) —> \GS)) 
and from an eigenstate at the center of the spectrum to 
the ground state (\C) — » \GS}), for both the Ising model 
and the LMG model. In Fig. |4] we show the final infi- 
delity for different system sizes TV as a function of the 
number of frequencies n/, at fixed total time T. In all 
cases considered, the infidelity decays exponentially with 
the rescaled number of frequencies np/B(N) = rif /TV a , 
showing a very similar behavior for both states in both 
models: indeed we have I ~ g(nf /N a ), g(x) being a scal- 
ing function of the form exp(— x 77 ), with 5 > 77 > 2 and 
1 < a < 1.5. The first interesting feature emerging from 
our analysis is that within each model the two transfor- 
mations I Ms) — >> \GS) and \C) —> \GS) approximately 
present the same complexity: a is only slightly larger 
for the I Ms) — >• \GS) conversion. This result can be ex- 
plained by the fact that the states \C) and \GS) are not 
trivially connected through the Hamiltonian, although 
they are both localized with respect to H[T] for T ^> 1. In 
practice also in this case the transformation is performed 
by first spreading the state onto the whole Hilbert space 
and then recombining the different amplitudes into the 
desired state. Such an operation requires approximately 
the same complexity as the state-to-state conversion be- 
tween the maximally spread state and the ground state, 
I Ms) — >• \GS). The second feature is instead emerging 
from the comparison between the two different models: 
the complexity scales approximately linearly with the size 
for both the LMG and the Ising model. 

Conjecture- The previous arguments suggest that the 
number of frequencies n f required in the optimal control 
pulse to achieve full control is not strongly dependent on 
the initial and final states but rather of the dimension of 
the manifold supporting the dynamics, namely D m (N). 
Here we conjecture that the number of independent basis 
functions rif needed to achieve a given level of fidelity, i.e. 
the complexity of the control problem, is rif ~ D m (N) 
up to subleading corrections. Since in both models con- 
sidered above the dynamics is supported by a manifold 
that scales linearly with the number of spins TV, in or- 
der to further support and verify our conjecture we will 
consider in addition a non integrable system, namely the 
Ising model in the presence of a longitudinal field J x 7^ 
in Eq. [TJ In this case, we expect that the optimization 
complexity should increase drastically since the effective 
manifold supporting the dynamics has dimension D m (N) 
now scaling exponentially with TV. We performed simu- 
lations for both cases, J x = (integrable system) and 
J x 7^ (non integrable system), analyzing the behaviour 
of the infidelity as a function of the number of frequen- 
cies for systems of different sizes. As before, in both 
cases the infidelity decays exponentially as a function of 
the number of frequencies, / ~ g(—rif/B(N)) (data not 
shown). In Fig. [5] we show the fitted decay rate values 
B(N) as a function of the size TV for J x ^ (full cir- 
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FIG. 5: Decay rate Basa function of the size TV for the Ising 
model with J x / (black circles) and J x = (red triangles). 

cles) and J x = (empty triangles). Despite the fact 
that, due to the exponentially growing Hilbert space, we 
are now restricted to small sizes 2 < TV < 8, as clearly 
shown by the fits, the rate B(N) for the integrable model 
(J x = 0) scales linearly with the size, while for the non 
integrable model (J x ^ 0) it scales as an exponential with 
TV, that is B(N) oc D m (N). This implies that in order to 
keep a constant infidelity g(—rif/B(N)) = const while 
increasing the system size TV, we need to scale the num- 
ber of frequencies like nf ~ B(N) ~ D m (N), in perfect 
agreement with the prediction of our conjecture. Inter- 
estingly, in the case J x = 0, optimal control complexity 
scales approximately linearly even though those results 
have been obtained simulating the Ising model in the 
full exponential Hilbert space without using the Wigner- 
Jordan transformation. Thus, optimal control complex- 
ity appears not to be influenced by the simulation details 
and thus CRAB control might be very effective in any in- 
tegrable system, even in cases where the exact solution 
is unknown. 

Conclusions- In this work we demonstrated for the 
first time the possibility of employing optimal control to 
reverse the dynamics of many-body quantum systems, 
effectively reducing the quantum entropy generated with 
strongly disordering processes. We showed that it might 
be possible to optimally reverse the system dynamics 
even in cases in which an exact reverse evolution cannot 
be realized (i.e. partially tunable Hamiltonians). Fur- 
thermore we demonstrated that the optimized reverse 
dynamics is extremely robust against external sources 
of noise. Finally, we proposed a conjecture establishing 
a relationship between optimization complexity and ef- 
fective dynamical dimension of quantum many-body sys- 
tems. We verified our conjecture in three different sys- 
tems, demonstrating that optimal control can work ef- 
fectively in exponentially large Hilbert spaces as long as 
the systems under consideration are exactly solvable (in- 
tegrable). Our results are relevant both for experiments, 
providing an alternative and more convenient implemen- 
tation of reverse dynamics, and to establish a connection 
among entropy, optimal control and integrability. We 
plan to extend this analysis to open system dynamics. 
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